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We present a broad range of measurements of the angular orientation 6a{t) of the large- 
scale circulation (LSC) of turbulent Rayleigh-Benard convection as a function of time. We 
used two cylindrical samples of different overall sizes, but each with its diameter nearly 
equal to its height. The fluid was water with a Prandtl number of 4.38. The time series 
6o(t) consisted of meanderings similar to a diffusive process, but in addition contained 
large and irregular spontaneous reorientation events through angles A9. We found that 
reorientations can occur by two distinct mechanisms. One consists of a rotation of the 
circulation plane without any major reduction of the circulation strength. The other 
involves a cessation of the circulation, followed by a restart in a randomly chosen new 
direction. Rotations occurred an order of magnitude more frequently than cessations. 
Rotations occurred with a monotonically decreasing probability distribution p{A0), i.e. 
there was no dominant value of A0 and small A0 were more common than large ones. For 
cessations p{A9) was uniform, suggesting that information of 6*0 (i) before the cessation 
is lost. Both rotations and cessations have Poissonian statistics in time, and can occur 
at any Oq. The average azimuthal rotation rate |^| increased as the circulation strength 
of the LSC decreased. Tilting the sample relative to gravity significantly reduced the 
frequency of occurrence of both rotations and cessations. 



1. Introduction 

The problem of Rayleigh-Benard convection (RBC) consists of a fluid sample heated 
from below [for example, see Siggia(1994), Kadanoff(2001), Ahlers et at (2002)]. The 
heat drives a convective flow and is thus transported out the top of the sample. In 
our case the sample is a cylindrical container filled with water. This system is defined 
by three parameters: the Rayleigh number R = agATL^/nv (a is the isobaric thermal 
expansion coefficient, g the acceleration of gravity, AT the applied temperature difference, 
L the height of the sample, k the thermal diffusivity, and v the kinematic viscosity), the 
Prandtl number a = ly/n, and the aspect ratio F = D/L (D is the diameter of the 
sample) . Convection happens as a result of the emission of volumes of hot fluid known as 
"plumes" from a bottom thermal boundary layer that rise due to a buoyant force, while 
cold plumes emitted from a top boundary layer sink. In the turbulent regime of F = 1 
samples that we study, these plumes drive a large-scale circulation (LSC) [Krishnamurty 
& Howard(1981), Sano et al. (1989), Castaing et al. (1989), Ciliberto et 0^.(1997), Qiu 
& Tong (2001a), Funfschilling & Ahlers (2004), Sun et al. (2005), Tsuji et al. (2005)], 
also known as the "mean wind" , which is oriented nearly vertically with up-flow and 
down-flow on opposite sides of the sample. 
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The LSC configuration does not have the rotational invariance of the cylindrical sam- 
ple, but the cylindrical symmetry implies that any azimuthal orientation of the LSC is 
an equally valid state for the system. In this paper we present extensive measurements of 
spontaneous angular changes of , i-e. of reorientations of the LSC. Such changes have 
been observed previously [Keller(1966), Welander (1967), Crcvcling et al. (1975), Gorman 
et aL(1984), Hansen et a/.(1992), Cioni et aL(1997), Nicmcla et al. (2001), Furukawa & 
Onuki (2002), Sreenivasan et al. (2002), Brown et al. (2005b), Sun et al. (2005), Xi et al. 
(2006)]. In one case, a rotation of the entire structure through an angle of ~ 1/2 rev- 
olution without a significant change in flow speed was clearly observed in an experiment 
using mercury as the fluid [Cioni et aL(1997)]. Another conceivable mechanism is cessa- 
tion, in which the LSC flow-speed vanishes, and then the flow restarts in a different direc- 
tion. Cessation was observed in numerical simulations [Hansen et oZ.(1992), Furukawa & 
Onuki (2002)], a dynamical-systems model [Araujo et al. (2005)], and a stochastic model 
[Benzi (2005)]. All of these cases are two-dimensional models where only cessations with 
= 1/2 revolution arc possible, although in principle the models could be extended 
to three dimensions. Cessation also occurs in convection loops (a thin circular vertically 
oriented loop filled with fluid heated in the lower and cooled in the upper half) where be- 
cause of the two-dimensional nature again only cessations with A6' = 1/2 revolution are 
possible [Keller(1966), Welander (1967), Crevehng et al. (1975), Gorman et a/.(1984)]. 

The experimental work by Niemela et al. (2001) and subsequent analysis by Sreeni- 
vasan et al. (2002) yielded statistics relating to reversals of the LSC, but could not 
determine A^ and was unable to distinguish between the rotation and cessation mech- 
anisms. Cessations were first documented in a laboratory sample of turbulent RBC by 
Brown et al. (2005b), and the present paper presents more extensive analysis and addi- 
tional data from that project. Recent experiments in F = 0.5 samples [Sun et al. (2005)] 
measured the orientation of the LSC, and contemporary experiments in F = 1 samples 
[Xi et al. (2006)] have produced results complementary to our own regarding the az- 
imuthal dynamics of the LSC. With these experiments, the azimuthal dynamics of the 
LSC are beginning to be understood. Spontaneous changes of the orientation of the LSC 
are not only interesting from the standpoint of fundamental physics, but are important in 
many geophysical applications. For instance, reversals occur in natural convection of the 
Earth's atmosphere [Doorn et al. (2000)] . Convection dynamics in the outer core of the 
Earth are responsible for changes in the orientation of Earth's magnetic field [Glatzmaier 
et al. (1999)]. 

The goal of the present work is to better understand the reorientations of the LSC, by 
both the rotation and cessation mechanisms. We first explain the experiment and how we 
determine in Sect. 2. In Sect. 3 we illustrate the existence and nature of both rotations 
and cessations. In Sect. 4 we present statistics relating to reorientations of the LSC, 
showing that successive rotations are independent of each other, and that they result in 
a wide range of A0, in which strict reversals (A0 ~ 1/2 rev.) are not especially common. 
In Sect. 5 we present statistics relating to cessations. These are found to be an order of 
magnitude more rare than rotations. We show that aftc;r the; LSC stops, it restarts with a 
random new orientation. In Sect. 6 we compare our results with those of Sreenivasan and 
coworkers. There we show that we can reproduce the statistics that they derived from 
their data only when we count an event each time when the orientation crosses a fixed 
angle. We refer to such events as "crossings" . Crossings include events caused by small- 
amplitude, high-frequency "jitter" near the crossing angle, and considering them resolves 
apparent inconsistencies between the two experiments. In Sect. 7 we present statistics of 
the azimuthal rotation rate over the long term. There we show that the angular distance 
travelled by scales as in a diffusive process, and that the absolute value of the rotation 
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rate |^| increases when the LSC amphtude decreases. In Sect. 8 we present results from 
tilting the sample relative to gravity. These data complement results already reported 
by Ahlers et al. (2005). Since naturally occurring convection systems are generally not 
cylindrically symmetric, it is important to study how rotations and cessations behave in 
less symmetric systems. As observed by others (see, for instance, Sun et al. (2005) and 
references therein), wo find that the tilt pushes the LSC into a preferred orientation along 
the direction of the tilt. In addition, we determined that both rotations and cessations 
are strongly suppressed by the tilt. Finally in Sect. 9, we comment on the contemporary 
experiments by Xi et al. (2006), who reported on the azimuthal dynamics of the LSC, 
and compare the results to ours. A brief summary is given in Sect. 10. 

Because the theoretical models [Araujo et al. (2005), Benzi (2005)] only predict two- 
dimensional reversals, and not the three-dimensional rotations and cessations that we 
observe, they cannot yet be compared in detail to the experimental data. 

2. The apparatus and experimental method 

The experiments were done with two cylindrical samples with aspect ratio F ~ 1 that 
arc the medium and large sample described in detail elsewhere [Brown et al. (2005a)]. 
Both had circular copper top and bottom plates with a plexiglas side wall that fit into a 
groove in each plate. There were no internal flanges, seams, sensors, or other structures 
that could interfere with the fluid flow. The medium sample had D = 24.81 cm and 
L = 24.76 cm, and the large sample had D = 49.67 cm and L = 50.61 cm. Each sample 
was filled with water and the average temperature between the bottom and top plates 
was kept at 40.0° C where a = 4.38. The two samples of different heights allowed us 
to cover a larger range of R at the same a and F, so the overall range studied was 
3 X 10^ ^ 10^^. Three rows of eight blind holes each, equally spaced azimuthally and 
lined up vertically with each other at heights 3L/4, L/2, and L/4, were drilled from the 
outside into the side walls of both samples. Thermistors were placed into them so as to 
be within d = 0.07 cm of the fluid surface. Earlier experiments were done with only the 
middle row of eight thermistors at height L/2, and presentations of data in this paper 
that mention only one row of measurements refer to the middle row, which was sampled 
in both the early and later experiments. Since the LSC carried warm (cold) fluid from the 
bottom (top) plate up (down) the side wall, these thermistors detected the location of 
the upflow (downflow) of the LSC by indicating a relatively high (low) temperature. No 
parts of the thermistors extended into the sample where they might have perturbed the 
flow structure of the fluid. The lead wires for these thermistors were wrapped around the 
insulating layers just outside the side wall to prevent the introduction of heat currents 
into the sides of the samples from these leads. The thermistors have a resolution of 10~^ 
degrees C. Both samples were carefully levelled to better than 0.001 rad, except for the 
experiments in which we deliberately tilted the samples. 

We presume that the side-wall thermistors measured the temperature of the thermal 
plumes and the accompanying LSC, just outside of the viscous boundary layer at the side 
wall. Here we address the issue of whether one can measure these temperatures through 
the boundary layer. The response time of the thermistors for thermal diffusion through 
the side wall is d^/Kgw = 0(1) s, where Kgw is the thermal diffusivity of the plexiglas 
wall. The response time for heat flow through the viscous boundary layer - assumed to 
have a width A = 0.464L x Re [Grossman & Lohse (2002) with flt parameter from 
Brown et al. (2006)] - for thermal diffusion is expected to be A^/k w 0.215L'^ / (kRc). For 
the large sample this time ranged from about 140s for AT ~ 1°C {Re — 2500) to about 
32s for AT ~ 20° C {Re — 11000). However, this calculation greatly overestimates the 
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Figure 1. (a): An example from the large sample of the temperatures at the horizontal mid- plane 
of the side wall as a function of the azimuthal angle 6 for R — 9.6 x 10^". Solid line: a fit of 
Ti — To + (5cos(j7r/4 — 9o),i — 0, . . . , 7 to the data. The fit yields the orientation 6o and an 
amplitude S that reflects the LSC. (b): The averaged normalized side- wall-temperature-profile 
(Tj — To)/ 5, sorted into bins with each bin covering a small range of ^ — do- Solid line: a cosine 
function. Vertical bar: typical sample standard deviation for each bin. 



real response time, since the data in this paper show that we could observe temperature 
changes that occurred over time scales as short as several seconds. Consistent with direct 
measurements of fluctuations in the viscous boundary layer [Qiu & Xia (1998)], the 
relatively fast thermal response time suggests that there is turbulent mixing in this 
boundary layer that enhances the heat transport to the side-wall thermistors. 

We made measurements with a sampling period 5t as short as 2.5 seconds, and fit the 
empirical function 

T, = To + (5cos(i7r/4- 6*0), i = 0, . . . , 7 , (2.1) 

separately at each time step, to the eight middle-row side-wall thermistor-temperature 
readings. An example of such a fit is shown in Fig. la. Deviations from a smooth profile 
are presumed to be due to the turbulent natiue of the system; for instance a hot plume 
passing by a thermistor will cause a higher than average temperature reading at that 
particular angular location. The fit parameter (5 is a measure of the amplitude of the LSC 
and 00 is the azimuthal orientation of the plane of the LSC circulation. As defined here, 
the orientation 9q is on the side of the sample where the LSC is warm and up-flowing 
and is measured relative to the location of thermometer zero, which was located on the 
east side of the sample. Typically the uncertainties for a single measurement were about 
13% for 5 and 0.02 rev. for Oq. Fitting to the cosine function does not yield a unique 
^0 because it is 27r periodic. The flnal value of was chosen as the one closest to of 
the previous timestep, thus allowing us to observe rotations of the LSC through larger 
angles, than if we had reduced the range to < 0o < ^tt. We calculated orientations 
dt and 9b and amplitudes 5t and 5b for the top and bottom rows by the same method 
as for the middle row. Throughout the paper, when we show data based on one row of 
thermistors it is for the middle row, while the top and bottom row measurements are 
only used when data is shown for all three levels. 

To test the validity of the sinusoidal fitting function, Fig. lb shows the side-wall 
temperature-profile normalized by the fit values. Each point is an average of the nor- 
malized termperature (Tj — Tq)/5 in a bin with a small range oi 6 — Oq. The standard 
deviation of the normalized temperature for each bin is about 0.30 and nearly indepen- 
dent oi 9 — 9q. It is shown as a vertical bar in the plot to indicate the typical size of 
temperature fluctuations. The data is in good agreement with a cosine function (solid 
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Figure 2. A time series of the orientation 6o and amplitude 6 showing rotations of the LSC. 
Data is from the medium sample at R — 1.1 x 10^°. Solid line: middle- row thermistors. Dashed 
line: top-row thermistors. Dotted line: bottom-row thermistors. 



line) without any additional fitting, indicating that Eq. 2.1 is a good function to represent 
the average temperature profile around the side wall relative to ^o- 

The side-wall thermistors are also used to obtain the plume turnover time. Autocor- 
relations of a single side- wall thermistor-temperature yielded peaks at times 0,T, 2T, 
etc., while cross-correlation functions of mid-plane thermistors on opposite sides of the 
sample were negative and yielded peaks at times T/2, 3T/2, etc. The peaks in the corre- 
lation functions indicate a periodicity in temperature fluctuations, i.e. plumes, circulating 
in the sample, hence we call T the plume turnover time. The technical details of this 
measurement will be reported by Brown et al. (2006). 



3. The nature of rotations and cessations 

We previously published time series of 6q and S for the samples with 8 side-wall ther- 
mistors [Brown et al. (2005b)]. Figure 2 shows a half- hour time series of and 6 from 
the medium sample with 24 side-wall thermistors at i? = 1.1 x 10^*^. The orientation 
and amplitudes are shown separately for the three rows of thermistors. The data contain 
a series of erratic rotations of the orientation of the LSC. Through these rotations, a 
reversal of the LSC direction (A^ ~ 1/2 rev.) is made over several hundred seconds, 
roughly during the time interval from 500 to 1200s. This is a slow reversal relative to 
the plume turnover time T, which is 49 seconds in this case. This type of reversal by 
rotation has been reported before, using temperature measurements at various azimuthal 
locations in the bottom plate of a convection cell to determine the LSC orientation [Cioni 
et a/.(1997)]. It is important to note that throughout these rotations, the temperature 
amplitude 5 remains non-zero. This implies that the LSC was circulating over this entire 
period, so these are rotations and not cessations. The values for the three rows arc seen 
to agree fairly well, which agrees with our assumptions about the vertical alignment of 
the LSC, but the top- and bottom-row temperatures generally have more variability than 
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Figure 3. A time series of the orientation 60 and of the amplitude 6 showing a cessation of the 
LSC, in which the amphtude drops to zero. Data is from the medium sample at i? = 1.1 x 10^°. 
Solid squares: middle-row thermistors. Open triangles: top-row thermistors. Open circles: bot- 
tom-row thermistors. 

the middle-row temperatures, and the amplitudes St and St drop on occasion without 
significantly affecting the middle-row signal. 

Figure 3 shows another time series of 9o and S for the three rows of thermistors. Here 
S decreased essentially to zero and then increased back up close to its average value. 
Wc interpret an amplitude drop to also indicate a velocity drop, since the temperature 
distribution - represented hy S - drives the LSC by buoyancy, and experiments have 
found a correlation between temperature and velocity [Niemela et al. (2001), Qiu et al. 
(2004)]. This means that the LSC gradually slowed to a stop, reversed direction, and 
gradually sped up again in another direction without significant rotation of the plane of 
circulation. This is clearly a cessation, but note that it is not strictly a reversal because 
AO < 1/2 rev. 

4. Reorientation statistics 

In this section we examine the statistics of reorientations, regardless of whether they 
occur by rotation or cessation. However, these results reflect primarily the properties of 
rotations because, as we shall show below in Figs. 7 and 14, cessations arc an order of 
magnitude more rare than rotations. We encountered a reorientation event roughly once 
per hour (see Fig. 7 below). The statistics of cessations will be discussed separately in 
Sect. 5. 

4.1. Definition of reorientations 
Because rotations do not necessarily have clear starting and ending points, and may 
have a wide range of sizes and speeds, they are difficult to define for the purpose of 
data processing and statistical analysis. Instead, we start by defining a reorientation as 
an event with a sufficiently large and quick change in the orientation of the LSC. More 
specifically, we required reorientations to satisfy two criteria, using only data from the 
middle row of thermistors since that was available from all of the experiments. These 
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Figure 4. (a): Probability distribution p(|A6|) of the angular change for reorientations. 
Open circles: experimental data. Solid line: power-law fit to the data, (b): The exponent 7 from 
the power-law fits of p(|A0|) as a function of R for the medium sample (solid diamonds) and 
large sample (open diamonds). Dashed line: the value of 7 from the fit to all data. 



criteria arc the same as those used by Brown et al. (2005b). First, the magnitude of 
the net angular change in orientation |A0| over a set of successive data points for 9q 
had to be greater than a chosen parameter J^Omin- Second, the magnitude of the net 
average azimuthal rotation rate \9\ = \/S.9//S.t\ over that set had to bo greater than a 
chosen parameter Omin- Here Ai is the duration of the reorientation. Usually multiple 
overlapping sets satisfied these requirements, so in those cases the set with the maximum 
local reorientation quality factor Qn = |A0|/(Ai)" was chosen as the reorientation. For 
< n < 1, Q„ represents a compromise between choosing the maximum angular change 
(Qo) or the maximum rotation rate (Qi). Any adjacent points to the chosen set were 
also included if the instantaneous rotation rate ^0 = SOo/St [56f) = 9Q{t + St) — Oo{t)] for 
the adjacent point was greater than 9min and of the same sign as for the reorientation. 
For the results presented in this paper, we used the parameters AOmin = 0.125 rev., 
9min = (0.1 rev. )/T, and n = 0.25; but since all three reorientation definition parameters 
are arbitrary, we did the analysis over the ranges 1/64 rev. < AOmin ^ 0.25 rev. , 
0.0125 rev. < 6minT < 0.2 rev. , and < n < 1 to confirm that the physical results were 
not sensitive to the choice of reorientation definition parameters. For AOmin and Omin, the 
smallest values used were on the order of the turbulent fluctuations in the data, so nearly 
every data point would be counted as a reorientation, and at the largest parameter values 
used too few events were counted to yield useful statistics. The qualitative conclusions 
of the analysis did not change with the different parameter values we tried. Some of the 
quantitative values found from fitting statistical distributions did change. Any variation 
of these fit values with the chosen parameters will be mentioned where it is appropriate. 

4.2. Results 

Probably the most important feature of reorientations is the angular c:liangc Ad. The 
probability distribution p(|A^|) is shown in Fig. 4a for all data, regardless of R. Here 
the data were sorted into bins that were evenly spaced on a logarithmic scale. The error 
bars represent the probable error of the mean, with the relative error of the mean taken 
to be the inverse square root of the number of reorientations in a bin. Fitting a power 
law p{\A0\) oc (|A6'|)'^ to the data yielded 7 = -3.77 ± 0.04. The fit was done by the 
maximum-likelihood method [sec, for instance, Bevington & Robinson (1992)] to avoid 
errors associated with the binning of the data. The same analysis was done also for re- 
orientations seperately at various R. These probability distributions were again fitted by 
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Figure 5. (a): The probability distribution p(|^|T) of the azimuthal rotation rate for reorien- 
tations at all R. Open circles: experimental data. Solid line: power-law fit to the data, (b): The 

exponent /i from the power-law fit of p(|^|T) as a function of R for the medium sample (solid 
diamonds) and large sample (open diamonds). Dashed line: the value of fj, from the fit to all 
data. 



power laws, with the resulting 7 values shown for each R in Fig. 4b. The exponent of 
the distribution was, within our resolution, independent of R. The power-law exponents 
found could vary by up to a factor of 2 for extreme values of the reorientation defini- 
tion parameters. However, for all parameter values tried, the distribution was consistent 
with a power law with a negative exponent that did not vary significantly with R, and 
there were never any peaks in the distribution. The major conclusions arc that there is 
a monotonically decreasing distribution of \A9\, so that smaller reorientations are much 
more common than larger ones, and that there is no characteristic reorientation size. 
This analysis shows that the strict reversal is not especially common among reorien- 
tation events. The probability distribution in Fig. 4 implies that only about 1% of the 
reorientations we counted have | = l/2±0.05 rev. This result is significant because an 
interpretation of previous experimental work had suggested that the events found were 
all reversals [Niemela et al. (2001)], and two-dimensional theoretical models can predict 
only reversals [Araujo et al. (2005), Benzi (2005)]. 

The average azimuthal rotation rate \9\ for reorientations was studied by a method 
analogous to that used for \A6\. All of the reorientations were sorted into bins according 
to \0\T, to make the rotation rate dimensionless so data at different R could be compared. 
The probability distribution p{\6\T) is plotted in Fig. 5a, with error bars ccjual to the 
probable error of the mean. Fitting a power law to the data in the range \0\T > 0.16 rev. 
yielded p{\9\T) oc {\9\T)'^, where = —2.61 ± 0.04. The same analysis was also done for 
individual values of _R, and /i is shown for each in Fig. 5b. The exponent /i can vary by 
about a factor of 2 with extreme values for the reorientation definition parameters, but the 
qualitative results are unchanged. The probability distribution p{\0\T) is described well 
by a power law with a largc^ nc^gativc exponent that is within our resolution independent of 
R, showing that slower reorientations are much more common than faster ones, and that 
there is no characteristic rotation rate for reorientations. We note that the normalization 
by T is not independent of R. In fact, T oc R~^-^ [Brown et al. (2006)] over most of our 
range of R. 

The probability distribution of the duration At of reorientations, which is not shown 
here, is sharply peaked, and the peak location coincides with A9rnin/9rnin- This is ap- 
proximately true for all of the values of the reorientation parameters we studied. Since 
this characteristic duration always depended on these artificial parameters, there seems 
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Figure 6. Solid lino: Probability distribution of the orientation of the mean wind p{9o) as a 
function of Oo. Open circles and dotted line; Probability distribution of the starting orientation 
p{0o,i) for reorientations. Solid circles and dashed line: Probability distribution of the ending 
orientation p{9oj) for reorientations. Data is for i? = 1.1 x 10^'^ in the medium sample. 



to be no characteristic physical time scale for the duration of reorientations that we could 
measure. Inspection of the data indicates that the duration of reorientations can be as 
short as O.IT (these appear to be cessations) and the longest durations seem to be about 
equal to the artificial value (1 rev.)/0„ij„. 

Since the method for defining and counting reorientations is non-traditional, the anal- 
ysis program was run on a simulated Brownian diffusive process for comparison. The 
orientation 9q was allowed to travel in one dimension in discrete steps at time intervals 
equal to the sampling interval St of experiments. The orientation change 69o for each 
time step was randomly generated with the distribution of Gaussian white noise. The 
simulated noise was made to have the same root-mcan-square step size as the diffusive 
region of the real data for i? = 1.1 x 10^° (given by ifVSi as defined in Sect. 7). For this 
simulated data there are about half as many reorientations as in the real data. While 
p(|A0|) looks similar in shape to that for the real data, it falls below the real data and 
an exponential distribution fits better than a power law to p(|A0|). This suggests that 
the LSC undergoes significantly more large reorientation events than a diffusive process. 

One notable point is that there is a preferred orientation 6m of the LSC, which is 
apparent in the probability distribution of the mean wind orientation p{9o) shown in 
Fig. 6 for i? = 1.1 X 10^° (solid line). In this case, we reduced the orientation to the 
range < 9q < 2tt. For other R, p{6o) is generally found to have a single broad peak at a 
9m that varies with R, but that is reproducible when experiments are done in the same 
apparatus at the same R. This distribution would ideally be uniform in an azimuthally 
symmetric system, but minor deviations from perfect rotational symetry such as a slight 
tilt of the sample relative to gravity, a slightly elliptical cross section of the side wall, or 
a coupling of the Earth's Coriolis force to the LSC [Brown & Ahlers (2006)] could cause 
a deviation from the uniform distribution. Also shown in the figure are the probability 
distributions of the starting orientations p{0o,i) (dotted line) and ending orientations 
p{9oj) (dashed line) of reorientations from the same data set. Both of these distributions 
fall reasonably close to p{9o). This shows that reorientations can occur at any orientation 
^0 of the mean wind. This conclusion differs from that of Niemela et al. (2001), whose 
interpretation of their data implied a bimodal distribution of the mean wind orientation. 
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Figure 7. The average rate of occurrence uir for reorientatfon as a functfon of the Rayleigh 
number R. Solid circles: medium sample. Open circles: large sample. 



with the orientation switching between two opposite orientations. This difference will be 
discussed further in Sect. 6. 

The average rate of occurrence Ur of reorientations is shown as a function of R in 
Fig. 7, with error bars equal to the probable error of the mean. This rate did not vary 
much with R for most of the range studied, but note that since Omin = (0-1 rev. )/T(R) , 
the minimum requirement is more stringent at larger R, so a definition for reorientations 
with 6min independent of R would result in the frequency of events increasing with R. 
The rate cOr depended strongly on the reorientation definition parameters: it decreased by 
about 2.1% with each increase of 1% in A6min, and it decreased by about 1.5% with each 
increase of 1% in Omin - When the new side wall with 24 thermistors with height L = 49.54 
cm was used as part of the large apparatus, uir increased by nearly a factor of 2. It was 
also found that p{9o) had a smaller peak and was closer to a uniform distribution (i.e. 
p{0 + 0) = 1) with the newer side wall, but other measured parameters such as Reynolds 
numbers did not change. The side walls were nominally identical, but the diameters had 
local variations of about 5 parts in 10,000. It is likely that the first side wall was less 
circular, resulting in a pressure that tended to force the LSC orientation to align with 
the long diameter of the side wall. This could explain why p{0o) had a stronger peak at 
the preferred orientation in the first side wall, and would suggest that this forcing of the 
LSC into a preferred orientation also suppreses reorientations. This suggests an incredible 
sensitivity of reorientations to the geometry of the system, and it is unclear why it would 
change the frequency of reorientations by such a large factor. This sensitivity does not 
seem to qualitatively affect other aspects of reorientations, such as p(|A6'|). While we can 
say that reorientations occur on the order of once per hour, it is difficult to draw more 
specific conclusions on this subject. 

We now consider the distribution of reorientations in time. Let r„ be the time intervals 
between the ith and (i + n)th reorientations. For n = 1, ti is simply the time interval 
between successive reorientations. This is similar to the definition used by Niemela et al. 
(2001), whose studies of the time intervals between events referred to by them as reversals 
of the LSC will be compared with the present work. However, their reversals were defined 
to occur at one instant, while in the present work reorientations are defined to have some 
duration. Thus we define ti as the time between, but not including the duration of, 
reorientations. The results are essentially the same when the analysis is done with the 
duration of reorientations included in n . 
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Figure 8. Solid circles: probability distribution p(ti/(ti)) of the time intervals n between suc- 
cessive reorientations for all R. (a): Over the entire range of n measured, (b): For n < (n). 
Solid lines: the function p(ti/(ti)) = exp(— ri/(Ti)), representing the Poisson distribution. Dot- 
ted line: fit of the power law p{ti/{ti)) oc {ti / {ti))~^ to the data, indicating a poor fit. 



All of the time intervals n were sorted into bins according to the value of t\/{t\), 
where (...) represents an average over a data set at a single value of R. The proba- 
bility distribution p(ti/(ti)) is shown in Fig. 8a over the full range of ti/{ti), and in 
Fig. 8b over the limited range n < (ri). The error bars indicate the probable error of 
the mean for each bin. The data are in good agreement with the exponential function 
p(ti/(ti)) = exp(— ti/(ti)), which represents the Poissonian distribution. Note that there 
are no adjustable parameters in this comparison. When reorientations follow Poissonian 
statistics, it means that successive reorientations occur independently of each other. 

The agreement of an exponential function with p{ti) for their reversals was found also 
by Sreenivasan et al. (2002), but only for large t\. They fit a power-law distribution 
p{ti) oc t^^ to their data for small time intervals, with ri < 1000s ~ SOT. This time 
interval roughly corresponds to t\ < (ti) in the present work. A fit of the function 
p(ri/(ri)) oc (ti/(ti))~^ to our data is shown as a dotted line in Fig. 8b for comparison; 
it is not a good representation of our data. For our results the exponential distribution 
p(ti/(ti)) = exp(— ti/(ti)) (solid line) continues to hold also at small ti. The cause of 
this difference between our results and those of Sreenivasan et al. (2002) will be discussed 
further in Sect. 6. 

As another test of the Poissonian nature of reorientations, an autocorrelation of suc- 
cessive time intervals is given by: 



{{Ti,k - {ri)n 

Here rijt is the kth time interval between successive reorientations when they are ar- 
ranged in order of occurrence. The plot of gT,r{n) is shown in Fig. 9 for i? = 1.1 x 10^°. 
According to the normalization, 5r,r(0) = 1 but for all other n > 0, gT,r{n) is scattered 
around zero. The autocorrelation function for a perfect Poisson process is a delta func- 
tion, while a finite sample size would result in some scatter around zero. Figure 9 shows 
good agreement with Poissonian statistics. 
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FiGLTRE 9. The autocorrelation of time intervals between reorientations gT,r{'n) for i? = 1.1 x 10"'^'^ 
in the medium sample. A Poissonian process should yield a delta function at n = and is 
consistent with the data. 



5. Cessation statistics 

During the entire investigation, spanning about one year of data acquisition in eacli 
of the two samples, we observed a total of nearly 1000 cessations. Of these 694 were for 
untilted samples and at the Prandtl number a = 4.38 under consideration in this paper. 
In addition, 52 cessation events were encountered for a = 4.38 in samples that were tilted 
relative to gravity at various angles. 



5.1. Definition of cessations 

To better distinguish between the rotation and the cessation mechanism, we now consider 
statistics for cessations only. In a previous paper [Brown et al. (2005b)] we identified ces- 
sations as a subset of reorientations by determining whether the amplitude 5 had dropped 
below a specified value during the reorientation. We found that cessations accounted for 
about 5% of reorientations. However, there also were events where the amplitude dropped 
without a significant change in orientation. These should be counted as cessations, but 
were not counted as reorientations. Thus we now redefine cessations and count them 
whenever 5 drops below a chosen minimum amplitude 5i. All of the adjacent points in 
the time scries arc counted as part of the cessation as long as 5 is below a chosen max- 
imum amplitude 5h > Si. The lower-amplitude threshold 6i was chosen as the largest 
value such that p(|A6'|) was uniform. It was found to depend on the side wall used, but 
for each side wall it was chosen as the largest value such that A0|) was uniform. For all 
side walls this limiting distribution of uniform p(| A^|) was reached for some small 5i and 
p(|A0|) remained uniform in the limit as Si was reduced to zero. For the medium sample, 
and for the large sample side wall with eight thermistors, we used the parameter value 
Si = 0.15(5). For the large sample side wall with 24 thermsitors, we used Si = 0.07(5). 
The upper-amplitude threshold Sh was chosen as the largest value such that S was equal 
to its limiting value for small S. For all samples this was 5^ = 0.5(5). These properties of 
cessations will be explained in more detail later in this section. Both parameters depend 
on the average amplitude (5) at each R, which increases with R. Because our definition 
of cessations depends on S while reorientations depended on Ad and 9, some events fall 
into both categories, possibly with different starting and ending times. Other events are 
unique to either reorientations or cessations. 
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Figure 10. The probability distribution p(|A^|) of the net angular change during cessations 
for all R. Solid Une: the uniform distribution. 



5.2. Results 

The probability distribution p(| A(?|) of the not angular change during cessations is shown 
in Fig. 10 for data at all R from both samples and including diffferent side walls with 
different Si. For cessations we calculated AO reduced to the range — tt < < tt by 
adding or subtracting multiples of 2tt, and then fiu'ther reduced it to \ A0\. The probability 
distribution is seen to be consistent with the uniform distribution. This agrees with our 
earlier results for p(| A0|) [Brown et al. (2005b)], where we counted cessations as a subset 
of reorientations. This current plot additionally covers the range < \A9\ < 1/8 rev., 
which could not have been counted using the reorientation algorithm. This p(|A0|) for 
cessations is very different from the distribution for reorientations, which was found to 
follow a power law with a large negative exponent. The uniform distribution of angular 
changes implies that after the LSC stops, it is equally likely to start up again at any new 
orientation, apparently losing its memory of its previous orientation. This interpretation 
is the reason for using the reduced range — tt < AO < tt, since there is no physical 
difference between choices of AO separated by 2tt for cessations, in contrast to the case 
for rotations where there is a continous variation in ^o- Again, reversals (i.e. AO =1/2 
rev.) are not especially common. 

Figure 11 shows the probability distributions p{9o,i) (solid squares) and p{Oqj) (open 
circles) of the orientations ^o,i and 6qj at the beginning and at the end of cessations 
for all data from the medium sample. The error bars represent the probable error of the 
mean of each bin. The distribution p{0o,i) at the beginning of cessations is peaked near 
a preferred orientation 6m, much like p{9o) (see Fig. 6). This distribution is non-zero 
for all 9q, so cessations can occur at any orientation. The distribution p{Ooj') at the 
end of cessations is consistent with a uniform distribution, thus it is consistent with our 
conclusion that the LSC restarts at a random orientation after a cessation. It is interesting 
to note that whatever inhomogeneity causes the maximum in p(0o) and p(0o,i) does not 
have any effect on p(9qj). This is consistent with the fact that cessations are relatively 
quick events, happening in about a turnover time, while the non-uniform p{6q) can be 
attributed to the net effect over a long period from a weak forcing [Brown & Ahlers 
(2006)]. 

One interesting quantity is the duration of cessations At. Figure 12 shows the proba- 
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Figure 11. The probability distributions p{6o,i) and p{9oj) ol the orientation at the beginning 
(solid connected circles, p{do,i)) and the end (open circles, p{Ooj)) of cessations. Dashed line: 
the uniform distribution. Data is for all R in the medium sample. 
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Figure 12. The probability distribution of the normalized duration of cessations p{At/T) for 

all R. 

bility distribution p{At/T) of At normalized by the turnover time. The figure shows a 
peak near At/T ~ 1.2. However, this peak location is dependent on the parameters we 
use in the defintion of cessations, since it represents the time for the amplitude to change 
from 5h, down to below Si, and back up again to Sh- 

To determine a more physically meaningful value relating to the duration of cessations, 
the duration of each cessation was calculated for several values of the; c;cssation cutoff 
parameter Sh- Each cessation was also divided into two time intervals, one before and 
one after the minimum amplitude was reached, so that the durations of the amplitude 
decrease and the amplitude increase could be determined separately. Figure 13a shows 
the average cessation half-duration (At/(2T)) for both the amplitude decrease (solid 
circles) and the amplitude increase (open circles) for different values of Sh/ {S). Cessations 
at all R were used so as to obtain a sufficiently large collection of events. The linear 
function {At/{2T)) = A{Sh - So)/{S) was fit to the data for Sh/{S) < 0.5. This yielded 
A = 1.61±0.04 for the amplitude decrease and A = 1.69±0.04 for the amplitude increase. 
It is interesting to note that on average the decay of the LSC during a cessation took 
just as long as the subsequent growth. The fit also gave So/{S) = 0.090 ± 0.005 for the 
amplitude drop and So/{S) = 0.099 ± 0.004 for the amplitude rise. The value of was 
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Figure 13. (a): The average normalized half-duration of cessations (At/(2T)) for differ- 
ent values of the cutoff parameter 5h / (S) for all 7?. Solid circles: duration of amplitude de- 
crease. Open circles: duration of amplitude increase. Solid line: fit of the linear function 
(At/(2T)) = A{5h — 6o)/{6) to the data with 6h/S < 0.5 for the amplitude decrease. Dashed 
line: fit of the same function to the amplitude increase, (b): The fit parameter A for several 
values of R in the medium sample (solid symbols) and the laxge sample (open symbols) for the 
amplitude decrease (circles) and increase (diamonds). Solid line: the average A over all R for 
the amplitude decrease. Dashed line: the average A over all R for the amplitude increase. 




Figure 14. The average rate of occurrence of cessations ujc versus R. Solid circles: medium 
sample. Open circles: large sample. Solid line: average frequency of cessations for all data. 



close to the average minimum amplitude for cessations {mm{S) / {S)) = 0.088 ±0.003, and 
probably represents a base level of temperature fluctuations without the LSC. 

Inverting the fitting function yields a characteristic rate of change of the amplitude 
\5\ = 2{Sh - 6o)/At = 1/A X {S)/T = 0.62 ± OM{S)/T for the amplitude decrease and 
0.59 ± 0.01 ((5) /T for the amplitude increase. Figure 13b shows how A varies with R for 
several data sets with at least 14 cessations each. This shows that the slope with the 
normalizations used above is independent of R, within the precision of the experiment. 
This justifies the use of data with all values of R in Fig. 13a. It is interesting that S has 
the same magnitude during both the decrease in amplitude and increase in amplitude, 
and remains constant for much of the duration of cessations. We had no reason to expect 
this, and this should put a significant restriction on any dynamical theories of cessations. 
The value of 6 also implies a physical time scale for cessations, so for a cessation with the 
amplitude starting and ending at 0.5(5), the duration is on average At = 2{6h — So)/\6\ « 
1.2T. 

The average rate of occurrence of cessations uic versus R is shown in Fig. 14 for data 
sets that are at least two days long. The error bars represent the probable error of the 
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Figure 15. The probability distribution of time intervals between cessations p(ti/(ti)) for all 
R. The entire range of the data is shown in (a), while a range restricted to n < (n) is shown 
in (b) for better resolution. Solid lines: p{ti/{ti}) = exp(— Ti/(ri)), representing the Poisson 
distribution. Dashed line: power law with an exponent of —1, showing a poor fit to the data for 
Ti < (n). 

mean for each data set. The data are consistc^iit with a rate lUc independent of R over 
the range studied, except for a decrease in uic for the smallest R. Over 226 days of total 
running time, the average rate was 1.49±0.08 day~^. It should also be noted that Wc does 
depend on the cessation definition parameter Si, which varied with the side wall used, 
although curiosly cOc did not change significantly with the side wall used. With each side 
wall considered separately, there is about a 1.7% increase in the cessation count with 
each 1% increase in Si/{S), although the value of 5i chosen is significant because it is the 
largest value that leads to a uniform p{A9). 

The probability distribution of time intervals between cessations p{ti/{ti)) for all 
data from the medium sample is shown in Fig. 15a. with error bars representing the 
probable error of the mean. The exponential function p{ti/{ti)) = exp(— ti/(ti)) agrees 
well with the data, indicating that cessations follow Poissonian statistics in time, as did 
reorientations. Figure 15b shows p(ti/(ti)) over the limited range ti < (ti), which shows 
that the data are consistent with the same exponential function (solid line) in this region. 
This distribution can not be fitted by a power law with exponent —1 (dashed line), as 
found by Sreenivasan et al. (2002) for their reversals. This difference between the two 
experiments will be examined further in the next Section. 

The autocorrelation of the time intervals between successive cessations gr,c{n), defined 
in a manner analogous to that given by Eq. 4.1, is plotted in Fig. 16 for i? = 1.1 x 
10^°. This data is from a run that contains 116 cessations over 54 days. Much like for 
reorientations, we find only scatter about gT,c{n) = for n > 0, providing more evidence 
that cessations have Poissonian statistics in time. 

6. Comparison of temporal statistics with earlier experiments 

An extensive study of LSC "reversals" was reported by Niemela et al. (2001). These 
authors used a F ~ 1 sample filled with helium gas at temperatures near 5 Kelvin. They 
reported data from a single pair of temperature sensors located at half-height in the 
fluid near the side wall. At that location they could determine the prevailing vertical 
component of the velocity of hot or cold temperature fluctuations, i.e. of "plumes" . It 
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Figure 16. The autocorrelation of time intervals between successive cessations gT,c{n) for 
R = 1.1 X 10^'^ in the medium sample. A Poissonian process should yield a delta function 
at n = and is consistent with the data. 



is generally held that the plume velocity is the same as that of the LSC. They defined 
reversals of the LSC as any switch between up-flow and down-flow of this local vertical 
velocity component. 

In a sequence of subsequent papers [Sreenivasan et al. (2002), Niemela & Sreeni- 
vasan(2002), Niemela et al. (2002), Sreenivasan et al. (2004), Hwa et aZ.(2005)] some 
of the members of the same research group carried out various statistical analyses of a 
particular data set taken at a Rayleigh number of 1.5 x 10^^, i.e. close to the highest 
Rayleigh number achieved in the present work . For this time series the Prandtl num- 
ber was 0.74 [Niemela & Sreenivasan(2003)], which differs somewhat from our a = 4.38. 
Although the value of a may have some influence for instance on the frequency of re- 
versals, it seems unlikely that the difference in a would qualitatively alter the physics 
of the reversals. The turnover time in our experiments is 49 s at i? = 1.1 x 10^*^, while 
Niemela et al. (2001) report a turnover time of about 30 s for i? = 1.5 x 10^°, so if the 
characteristic time scales for reversals are proportional to the turnover time, then we 
should expect them to be a little shorter for Niemela et al. (2001). 

The major difference between the results of their work and of ours is the distribution 
of p(ri/(ri)) for the time intervals ti between successive events. Sreenivasan et al. (2002) 
found that a power law fit the distribution p(ti/(ti)) for small ti and an exponential 
function fit the tail, while we found an exponential function to fit the distribution over 
the entire range for both reorientations and cessations (see Figs. 8 and 15). In a related 
ma.ttc^r there is a large difference between the frequency of events determined in the 
two experiments: Sreenivasan et al. (2002) reported about 420 events per day while 
we find only about 50 reorientations per day - an order of magnitude less. If their 
reversals correspond to cessations, then the difference is 2 orders of magnitude; or if the 
correspondence is to reorientations, for instance with A0 = 1/2 ± 0.05 rev., then the 
difference is 3 orders of magnitude. 

These differences can be understood by reinterpreting some results. Niemela et al. 
(2001) reported a nearly bimodal velocity distribution, which led them to the then rea- 
sonable belief that the large-scale circulation-plane had a more or less fixed azimuthal 
orientation relative to their sensors, and that the LSC was switching between two oppos- 
ing directions in that plane. Our experiments with eight side-wall thermistors provided 
new information in the azimuthal dimension, showing that the LSC in our system sam- 
ples all azimuthal orientations, and that reorientations and cessations can occur at any 
orientation. While it is possible that an asymmetry of the apparatus used by Niemela 
et al. (2001) could have resulted in two opposing preferred orientations, it seems unlikely 
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Figure 17. p(ti/(ti)) for crossings for R = 1.1 x 10^° in the medium sample, defined as occurring 
when the orientation of the LSC crosses either angle orthogonal to the preferred orientation drn- 
Solid line: a fit of the power law p{ti/{ti)) oc ti/(ti)~^ to the data for n < (n). Dashed line: 
a fit of an exponential function to the data for n > (ri). 

that this could also account for the much greater frequency of events observed by them, 
especially since our experiments with an asymmetry introduced by tilting the samples 
(see Sect. 8) show a strongly reduced frequency for both reorientations and cessations, 
but otherwise no qualitative difference. 

Here we present an alternative explanation for the results of Niemela et al. (2001). 
We consider the velocity in the plume layer, a thin region between the bulk and bound- 
ary layer near the side wall. We hypothesize that this velocity, when averaged over an 
intermediate time interval (say of order T), has a vertical component that is nearly a 
step function of 6, with speed +v when 9 is within ±7r/2 rad of 6q and —v for all other 
angles. To be clear, this time average should be over a period long enough to identify the 
LSC, but short enough to not contain any reorientations. The one-turnover-time average 
used in the experiment of Niemela et al. (2001) satisfies these criteria. A justification for 
this step-function velocity-distribution, independent of its ability to explain the results of 
Niemela et al. (2001), is given in the Appendix. As a consequence of random azimuthal 
meandering of the LSC this distribution would produce the; biniodal velocity distribution 
reported by Niemela et al. (2001). It also leads to an alternative interpretation of the 
events they counted: if reversals are counted whenever the local velocity switches direc- 
tions at a sensor as seen by Niemela et al. (2001), then this is equivalent to counting 
events whenever the LSC orientation crosses a critical angle orthogonal to the azimuthal 
sensor orientation so that the interface between the up- and down-flow moves past the 
sensors. These events usually arc not reversals with A9 = 1/2 rev., as they would include 
events covering a continuum of A9. Indeed for many of the events the LSC orientation 
would change only slightly, but just enough to cross the critical angle. We can test this 
hypothesis by reanalyzing our own data to study this new type of event which we shall 
call a crossing. In our system we shall count it whenever the LSC orientation 9o{t) crosses 
either of the two angles orthogonal to the preferred orientation, i. e. 6mi'n'/2 rad. This 
choice of the crossing angles is arbitrary as long as they arc separated by tt, and if we use 
for instance 9m and 0^ + 77 rad then our statistical results are qualitatively unchanged. 
If our hypothesis about the velocity distribution is correct, then the analysis of crossings 
from our data should yield the same statistical results as the events reported by Niemela 
et al. (2001) and Sreenivasan et al. (2002). 

The probability distribution of the time intervals p{ti/{ti)) for crossings is shown in 
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Figure 18. p(ri/(ri)) for crossings at 7? = 1.1 x 10^'^ in the medium sample, but now there is 
an extra buffer region of width 7r/4 rad around each crossing angle that the orientation must 
also exit before a crossing is counted. Solid line: a power law fit of p(ti/(ti)) oc (t"i/{ti))~^ to 
the data for ti/(ti) < 1. Dashed line: the exponential function of p(ti/(ti)) = exp(— ri/(ri)), 
representing the Poisson distribution. 

Fig. 17 at i? = 1.1 X 10^°, with error bars equal to the probable error of the mean as 
before. For one data set we found 2032 crossings over 11.8 days (for comparison wc found 
555 reorientations and 21 cessations in the same data set), corresponding to (ri) = 502 
s or about 172 crossings per day. This is only a factor of 2.4 smaller than the frequency 
of events reported by Sreenivasan et al. (2002), and thus plausibly consistent, especially 
since their experiment has a somewhat shorter turnover time. 

The figme shows a fit of an exponential function p{t{) oc exp[— ti/(/i(ti))]) to the data 
for Ti > (ri) (dashed line), which fits well in that range of ti but falls well below the 
data for smaller n. We obtain /i(ti) = (3.3 ± 0.1)(ti) = 1600 s from this fit. In the other 
limit, a power law p(ti/(ti)) cx {ti/{ti))~^ was fit to the data for n < (n) (solid hue), 
which fits well in that range of ri but is much higher than the data for larger ti . These 
results of both the frequency of events and of p{ti) for crossings are in good agreement 
with the results of Sreenivasan et al. (2002), which implies that the two experiments are 
consistent, but that crossings rather than rotations or cessations are comparable to the 
events reported in the earlier work. 

To understand why the power-law distributions occur, we again reanalyze the present 
data with a modified definition for crossings. It was suggested by Sreenivasan et al. 
(2002) that the LSC orientation undergoes some "azimuthal drift (or jitter)" . We want to 
eliminate the counting of events in which the orientation jitters back and forth around the 
crossing angles, and only count events where there is a significant change in orientation 
of the LSC. To this avail we define buffered crossings as occurring when 9o{t) not only 
crosses Om^T^/^ rad, but also exits a buffer region of width 7r/4 rad centered about either 
crossing angle. The probability distribution p{ti/ (ti)) for buffered crossings is shown in 
Fig. 18 for R= 1.1 X 10^°. An exponential function p{ti) = exp(— ti/(ti)) is consistent 
with the data (dashed line), showing that these buffered crossings are Poissonian. A 
power law p{ti) oc is also shown (solid line) for comparison. Now it is apparent that 
jitter, or small orientation changes of the LSC around the crossing angles, was responsible 
for the power-law distribution of p(ri) for small ri. There are only 520 buffered crossings 
from the same data set that had 2032 crossings, so the jitter is of course also the reason 
for the large number of events. 

We can continue the analysis of the temporal statistics for both reorientations and 
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Figure 19. Number of non-empty bins Nr from the binning of events on the time axis with 
bin width r for Ji = 1.1 x 10^" in the medium sample. Solid circles: reorientations. Open circles: 
crossings. Dashed lino: prediction for a Poisson process. Solid line: power law fit of Nr = r~^ to 
crossings for r < 0.3, yielding the fractal dimension D. Dotted line: space-filling limit Nr oc r~^ 
for crossings. 



crossings in the spirit of Nicmcla et al. (2001) and Srccnivasan et al. (2002). One aspect 
of that analysis was the fractal dimension of the time series of events. We took a time 
series with reorientations, divided up the time axis into bins of width r, and then placed 
each reorientation into the appropriate bin based on the time when the event occurred. 
We then counted the number of non-empty bins Nr- The fractal dimension D is defined 
by the equation Nr = r~^ in the limit of small r. Figure 19 shows Nr for various bin 
widths r for reorientations at all R (solid circles) and for crossings at i? = 1.1 x lO"'^" 
(open circles). For r :» (ti), the data always follows Nr oc r~^, which is the case where 
the bin width is so large that every bin contains at least one reorientation. This is called 
the space-filling limit. For r <C (ti), Nr reaches a constant (D = 0) for reorientations. 
This is the limit where the bin widths are so small that every reorientation falls into 
its own bin. Also shown in Fig. 19 is the analytical result for a Poissonian distribution: 
Nr — (No{Ti)/r) X [1 — exp(— r/(Ti))], where A^o is the total number of events. The 
theory agrees perfectly with our data for reorientations, so it provides further evidence 
that reorientations are Poissonian. However, the data for crossings can be fitted by a 
power law Nr cx r^^ for r < (ri), and yield D = 0.21. This is reasonably consistent 
with the results of Niemela et al. (2001), which gave a range of values 0.2 < D < 0.5 
for different R. The non-zero value of D is characteristic of the power-law distribution of 
time intervals; no matter how small r is, there are a finite number of time intervals that 
are smaller than r, so Nr will not reach the total number of events. 

When we calculate the fractal dimension for buffered crossings we get D = 0. D is 
indistinguishable from zero for buffer widths greater than Stt/IG rad., while D gradually 
increases for smaller buffer widths up to £> = 0.21 for no buffer. This suggests that once 
the orientation is 37r/32 rad. from the crossing angle it is far enough away that the jitter 
of 00 no longer affects the event statistics. 

The power-spectral density of the detrended timing of events was reported by Sreeni- 
vasan et al. (2002). Let the time of the k^^ reorientation tk from the beginning of a data 
set be defined as tk — {Ti)[k + A{k)]. Here A{k) is the deviation from a linear trend of 
events in time, where the linear trend is (ti)A; as defined by Srccnivasan et al. (2002). The 
power-spectral density E{u)) of A(fc) is given by E{uj) = l/(27r)| X;fc A(i)e-*'^'=|2. E{u}) 
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Figure 20. The power spectral density E{ijj) of the detrended time A(fc) of the fc*'' event for 
crossings (upper data) and reorientations (lower data) at R= 1.1 x 10^" in the medium sample. 
E(ijj) for rcorcintations is shifted down by a factor of 100. The solid straight lines are power 
laws with exponents of —2. 
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Figure 21. Moments (r') of the generalized time intervals between reorientations versus nth 
interval for q = 1,2,3,4,5,6 at 7i = 1.1 x 10^" in the medium sample. Ascending data sets 
correspond to increasing q. Solid lines: power law fits for each q of (t^) oc n'" for n > 100. 
Dotted lines: power law fits for each q of (r^) oc n*" for n < 5. 



is plotted in Fig. 20 for both crossings (upper data) and reorientations (lower data). In 
both cases, E{lo) is in good agreement with an w"^ roll off, consistent with a Brownian 
process, and consistent with the results of Srccnivasan et al. (2002). 

The moments of the generalized time intervals between events are given by (r^) = 
{{tk+n — tkY)k- The first six moments 1 < g < 6 are shown in Fig. 21 for reorientations 
for i? = 1.1 X lO"'^" in the medium sample. These data can be fitted by power laws 
(r^) (X n'» in limited ranges of n with different C for the large and small limits of n. 
The C values are shown as a function of q in Fig. 22. The small-n exponents (solid 
circles) were obtained from fits in the range n < 5, while the largc-n exponents (open 
circles) were from the range n > 100. For Poissonian data, we would expect ( = q for 
all n. This is in agreement with the large-n fits. The small-n data deviate significantly 
from the Poissonian prediction for large q. Srccnivasan et al. (2002) reported a similar 
deviation for their events from the Poissonian prediction. Together with the exponential 
cutoff of their power-law result for p{ti), they explained this as a finite-size effect. That 



22 



Erir Drown aitd Gii enter Aiders 




Figure 22. Exponents ^ vs. q from the fits (r^) oc n*" for n < 5 (solid circles) and for n > 100 
(open circles) for ii = 1.1 x 10^° in the medium sample. The solid line shows the prediction for 
Poissonian data. 
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Figure 23. Exponents ( vs. g from the fits (r^) cx n"- for crossings at _R = 1.1 x 10^° in the 
medium sample. Solid circles: for n < 5. Open circles: for n > 100. Dashed line: linear fit to 
data for n > 100. Solid line: quadratic fit to data for n < 5. 



explanation will not work for us because we found p(ti) to agree with an exponential 
distribution over the entire range of n. Inspection of p{ti) in Fig. 8 shows more large 
time intervals than would be expected for a Poisson distribution, although only by a few 
events. This might account for the large higher-order moments that we found, and thus 
for the offending values of (. The question remains as to whether reorientations are a 
Poissonian process and we recorded a lot of large time intervals by chance, or whether it 
indicates non-Poissonian physics at large time intervals. There are very few data points 
with these large time intervals, so we can not distinguish between these two possibilities. 

Again, to better compare to the results of Srccnivasan et al. (2002), wc calculated the 
moments for crossings (not shown) and C,{q) as before. Figure 23 shows the exponents 
C,{q) for both the large-n and small-n limits for crossings, but iov q < 1. For n > 100, 
^ cx g as would be expected for Poissonian data. For n < 5, the quadratic function 
Q = aq + bq^ was fit to the data and yielded a = 1.69 and b = —0.692. Similarly, a fit of 
a quadratic function for ({q) to their data was reported by Sreenivasan et al. (2002) and 
attributed to an underlying log-normal distribution. 

Another result from Sreenivasan et al. (2002) was that p(r„) for an intermediate num- 
ber of gaps, with 5 < n < 100, agreed with a bilognormal distribution. For crossings from 
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Figure 24. The probability distribution of the time intervals between 20 crossings pij-zo) for 
R= 1.1 X 10^° in the medium sample. Solid line: fit of a bilognormal distribution to the data. 



the present data, p([logr2o]^) is plotted in Fig. 24. A bilognormal function given by 



P([logr2o]^) oc exp 



[(logT2o)'-m]^ 



gave a good fit to our data. Again, this shows good agreement between our crossings and 
the reversals of Sreenivasan et al. (2002). 

In this section we showed that we can reproduce the major statistical results of Niemela 
et al. (2001) and Sreenivasan et al. (2002), but only when wc consider crossings. The 
results from our experiment differ from those of the earlier work when we consider only 
events corresponding to our definition of reorientations or cessations. The events reported 
by Nicmcla et al. (2001) were assumed by them to be reversals of the LSC because the 
local flow velocity reversed directions. Our results, based on the angular information 
about the LSC orientation afforded by the use of eight side-wall thermistors, have shown 
that they were likely counting crossings in which the orientation of the LSC crossed a 
fixed line, which often counts small jitters of the LSC orientation as events. Since we 
are interested in events that involve larger changes of the flow field, we have focused on 
studying reorientations consisting of rotations and cessations. 



7. Azimuthal rotation rate 

The instantaneous azimuthal rotation rate 9q is an important parameter relating to 
reorientations, and some of its statistical properties are covered in this section. The 
probability distribution pd^ol) of the absolute value of 6^ is shown in Fig. 25a for R = 
1.1 X 10^°. The function p(|^o|) = Sq/{\ + ax^) is fit to the data. The exponent £, 
representing the power-law dependence of the tail of p(|^o|)) is shown for various R in 
Fig. 25b. On average it is about 3.6 and it does not appear to vary much with R. 

The half width at half maximum of p(|^o|) is given by u; = a~^/^, but before we 
can normalize it to get a Reynolds number, we must take into consideration that av- 
eraged values of 9q depend on the experimental sampling interval 5t because the rota- 
tion is dominated by fluctuations. To see this dependence, the root-mean-square ro- 
tation rate 6™^ was computed for integer multiples of the sampling period n5t as 
= {[e^lt + nSt) - en{t)Y)/{n5t). This is shown vs. n5t as open circles in Fig. 26a 
for i? = 1.1 X 10^" for a carefully leveled sample (/3 = 0, see Sect. 8 below). The function 
Qrms _ ^p . (^Ji)-i/2 jg shown to fit the data over a large range of n5t (solid line). The 
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Figure 25. (a): The probability distribution of the absolute value of the instantaneous az- 
imuthal rotation rate p(|^o|) for i? = 1.1 x 10^" in the medium sample. Solid line: a fit of 
p(|^o|) = 'S'o/(l + ax^) to the data, (b): The exponent e representing the power-law dependence 
of the tail of the distribution vs. R for the medium sample (solid circles) and the large sample 
(open circles). 




n6t (sec.) 



Figure 26. (a): Root-mcan-squarc rotation rate 9™" as a function of the time interval nSt for 
R = 1.1 X 10^" in the medium sample. Open circles: data for a level sample (/3 = 0) with a 
fit of oc {n5t)~^^^ (solid line). Open squares: data for a tilted sample (see Sect. 8) with 

13 = -0.026 rad with a fit of oc {n5t)~'^ (dashed lino), (b): Reynolds numbers for the 

rotation rate vs. R. Circles: _Rf for 9™" , with a power-law fit (solid line) yielding the exponent 
X = 0.278. Diamonds: R"^ for the width w of p(|^o|), with a power-law fit (dashed line) yielding 
the exponent 0.245. Solid symbols: medium sample. Open symbols: large sample. 



{n5t) ^1"^ dependence indicates a diffusive process, i.e. Gaussian distributed fluctuations 
dominate Q™^ in this range of large n^t. 

We wish to define a Reynolds number ii^ = Li^O™^ /v, but we must choose a 6™^ 

corresponding to a standard time interval, for instance the viscous diffusion time tv = 
iP' jv^ so that data from different samples can be compared. This choice yields ^J""** = 

• and thus ijf = ^Lj^. This Reynolds number is plotted in Fig. 26b vs. R 
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Figure 27. (a): The average rotation rate as a function of amplitude for i? = 1.1 x 10^" in the 
medium sample with a fit of the power law (l^ol/(i^ol)) = rn{5/{S))^^ to the data (solid line), 
(b): The coefficient m (triangles) and exponent n (diamonds) vs. R from the power law fit for 
the medium sample (solid symbols) and the large sample (open symbols). 

(circles) and a power law i?f = cR^ is fit to the data (solid line) to obtain c = 0.0211 
and X = 0-278. 

We now return to the width w of p(|^o|)- In this case w was calculated only for the 
experimental sampling interval ^t, but we want to scale it up to assuming the same 
[nbt)'^!"^ scaling. Thus wc multiply by \/&tJtZ and define R^ = l^wjv x ^Jbtjtl = 
wLsJdt/v. This is also shown in Fig. 26b (diamonds). This scaling again allows us to 
compare Reynolds numbers for different experiments with different L. Indeed we see that 
the large and medium samples yield results that follow the same power law. However, 
i?™ does not indicate the actual rotation rate over the time scale t^. The latter is better 
represented by 6™^. The power law R^ = fR^ was fit to the data (dashed line) to obtain 
/ = 0.025 and g = 0.245. Wc sec that g is quite close to x- 

Next we consider the relationship between the magnitude of the rotation rate 16*01 and 
the amplitude 5. All of the |^ol for all R were normalized by the time average (|^ol) for 
their R and sorted into bins according to the value of the similarly normalized amplitude 
6/ {6) at the same time step. The average value of the rotation rate as a function of the 
amplitude is plotted for each bin as (|^o|/(|^o|)) vs. 5/ (5) in Fig. 27a for i? = 1.1 x 10^°. 
The bars represent the sample standard deviation of \9q\ for each bin, indicating the 
typical range of 6o. A fit of the power law(|0o|/(|^o|)) = m{5 / {5))~^ to the data yielded 
m = 0.908 ± 0.009 and tt = 1.16 ± 0.06. The same fit was also done separately for several 
different values of R and the values of m (triangles) and tt (diamonds) are shown in 
Fig. 27b. They appear to be independent of R. This analysis shows a nearly inverse 
relationship between the amplitude and the rotation rate. While a negative correlation 
between the two parameters was already evident during cessations, this shows a more 
general relationship. This is also consistent with a correlation function between 9q and 
5 published before [Brown et al. (2005b)], which showed a strong negative correlation 
between the two parameters. 

8. Tilting the sample 

The effects of tilting the sample relative to gravity were investigated before [Ahlers 
et al. (2005)]. Here we give some additional data relevant to the influence on reorienta- 
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Figure 28. Solid circles: The frequency of reorientations uir (left ordinate) versus tilt angle (5 
form the medium sample at i? = 1.1 x 10^°. Solid line: Gaussian fit to medium sample data. 
Open diamonds: the frequency of cessations u!c (right ordinate) vs. /? in the medium sample at 
i^= 1.1 X 10^°. 



tions. Tilting the sample by an angle f3 breaks the azimuthal symmetry and encourages 
the fluid in the thermal boundary layers to flow up the bottom plate and down the 
top plate in the direction 6^ of the steepest slope. Thus, 9j3 is also the preferred ori- 
entation of the LSC in a tilted sample when there are no stronger asymmetries. First 
we consider the effect of tilting the sample on the rate of occurrence of reorientations. 
Since T varies somewhat with /3 [Ahlers et al. (2005), Brown et al. (2006)], we used 
the turnover time for /3 = to determine Omin for counting reorientations for all of the 
tilting experiments. Figure 28 shows the average rate of occurrence of reorientations LOr 
(solid circles) versus tilt angle /3 at i? = 1.1 x 10^° in the medium sample. The error bars 
represent the probable error of the mean. The analogous data for R = 9.4 x 10^*^ and 
the large sample was shown by Ahlers et al. (2005). For both values of R the reorienta- 
tions are suppressed significantly even for very small tilt angles. The Gaussian function 
u>r = (27rcr5)~^''^ exp[— (/? — 6(3)^ /{2a^)], chosen empirically, was fitted to the data and 
gave a standard deviation of = 0.0160 ± 0.0006 rad, indicating how small of a tilt 
is required to significantly reduce the number of reorientations. For R = 9.4 x 10^°, the 
standard deviation was a^; = 0.0178 ±0.0011 [Ahlers et al. (2005)], indicating that there 
is no significant i?,-dependenee for the suppression of reorientations by tilt. The fit of 
oJrif}) is symmetric around a center offset of SP = 0.0022 ± 0.0006 rad, instead of (5/3 = 
as might have been expected. We will show elsewhere that this offset is consistent with 
the infiuence of the Earth's Coriolis force [Brown & Ahlers (2006)]. 

Figure 28 also shows the average frequency of cessations uic (open diamonds) versus 
/3 for three tilt angles where we have more than 10 days of data, all at i? = 1.1 x 10^°. 
The error bars represent the probable error of the mean. We also took data for 11 days 
total at larger tilt angles, with 0.017 rad < |/3| < 0.21 rad, and found no cessations in 
that period. Cessations are significantly suppressed in a tilted sample, about as much as 
reorientations, but because of the scarcity of events we cannot say more on this subject. 

The characteristic azimuthal rotation rates were also found to decrease in a tilted 
sample. For instance, the instantaneous rotation rate |^o| was reduced to about 50% of 
its level-sample value at /3 = —0.21 rad. This is not unexpected given the increase of the 
amplitude 6 by about 75% over its level-sample value for (3 = —0.21 rad [Ahlers et al. 
(2005)], and the inverse relationship between the two parameters shown in Fig. 27. 

For more than a minimum tilt angle, the root-mean-square rotation rate 0™* no 
longer scales as {nSt)~^^^ for large nSt. Rather, for |/3| > 0.026, it scales as {n5t)~^ for 
large nSt. This is shown in Fig. 26 for ii = 1.1 x 10^° and /3 = —0.026 (open squares) 
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Figure 29. The saturation value 59max as a function of the width ao of p{9o), for tilt angles 

> 0.026 in the medium sample at R = 1.1 x 10^". The solid line is the theoretical result 
S9max = \/2(7e that should pertain if SOmax is limited by the Gaussian distribution of p(^o)- 

with a fit of = S9max/{nSt) to the data, where d9max is a suggestively named fitting 
parameter. The {nSt)~^ scaling indicates that (Jf?™" saturates at a maximum value S6max 
at n6t « 300 s. Presumably this occurs because 9q is locked into a small range for large /3, 
which would suggest that 60max is related to p{9o). For tilt angles |/3| ^ 0.0087, p{9o) is fit 
well by a Gaussian distribution with width [Ahlers et al. (2005)]. For an ideal Gaussian 
distribution, the mean-square distance between two points x and y in the distribution is 
given by 

dx / dy{x - yf{2Tral)-' eM-^Vi^al)] eM-yyi2al)] = 2al 

-OO J — C30 

Thus if the locking of into a small range for large /3 is responsible for the saturation 
of (56',™% wc should have 59 max = v^ce- Figure 29 shows 59„iax vs. tJo- The data point 
on the upper right is for /3 = —0.026 and the point on the lower left is for /3 = —0.21. 
Also shown in the plot is a solid line for the theoretical prediction 56 max = V^ae. There 
is excellent agreement between the two, indicating that the large n5t dynamics of 9q for 

> 0.026 are dominated by the tilt locking the orientation into a small range, which 
overwhelms the diffusive dynamics that dominate the large n5t range for /? = 0. Even 
for very small /3 = —0.0044 we see deviations from the (nSt)'^^"^ scaling, indicating that 
these diffusive dynamics only dominate in very symmetric systems. 

It is interesting to consider why reorientations are suppressed by tilting the sample. If 
tilting simply added a slight deterministic rotation rate due to the added buoyancy of 
the boundary layers, then we would expect it to suppress reorientations with |^| smaller 
than the deterministic term. We can use a highly simplified model to estimate the order 
of magnitude of this effect. This model is based on one used successfully to estimate 
the effect of tilting on the Reynolds number [Chilla et aL (2004b), Ahlers et al. (2005)]. 
Tilting the sample by a small angle (3 relative to gravity results in a buoyancy force per 
unit area in the thermal boundary layers approximately given by plgl3aAT/2 parallel to 
the plates, where / is the boundary-layer thickness. This forcing mainly contributes to 
enhancing the LSC since it is usually aligned with the slope of the plates [Ahlers et al. 
(2005)], but when the LSC is not aligned with this slope, some fraction of this forcing will 
push the LSC into an azimuthal rotation towards the orientation 9/3. We would expect 
this fraction to be proportional to sin(6'^ — ^o)- This buoyant forcing is opposed by the 
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Figure 30. The average azimuthal rotation rate (^o(^o, Tmax)} with a delay time Tmax after the 
LSC has reached the orientation Oo. Medium sample, 7? = 1.1 x 10^° and /? = —0.0044. SoUd 
line: fit of {9o{0o,Tmax)) = Asm{df} — 9o) to the data. 



viscous shear stress from the azimuthal motion across the boundary layer. This opposing 
force can be approximated by pvue/l ~ pvLOo/ {21). Equating these terms, substituting 
I = L/{2M) {Af is the Nusselt number), and using the definitions of R and a yields 

. l3Rusm{00-eo) 
■ (^-^^ 

We carried out experiments at i? = 1.1 x 10^° and a = 4.38 {Af = 133, independent 
of (3 to better than 1% [Ahlers et al. (2005)]) at various tilt angles. Substituting these 
values into Eq. 8.1 we obtain Oq = 0.40/3sin(^^ — ^q) rad/s. By expressing 9q as a function 
of ^0 we implicitly assumed that inertia is negligible, i.e. = 0. This model also ignores 
possible buoyant forcing in the bulk due to the non-uniform temperature distribution 
there, and additional viscous stress in the viscous boundary layers as opposed to just 
the thermal boundary layers which were chosen for the convenience of balancing stresses, 
as well as any deformation of the LSC due to uneven distribution of these forces. This 
model is meant to predict a shape for the distribution of 0(){9q) due to buoyant forces, as 
well as the order of magnitude of the effect. Neither of these predictions are expected to 
change significantly by accounting for the aforementioned imperfections of the model. 

We can find the magnitude of this effect experimentally by averaging the instantaneous 
rotation rate at different orientations ^^o- Even though the turbulent meanderings of 
the LSC orientation tend to be much larger than the deterministic part, by averaging over 
many data points at the same value of we arc able to resolve the deterministic effects. 
The data were binned according to ^o, and the average azimuthal rotation rate {6q{6o)) 
was calculated for each bin. However, a similar calculation for {0o{9o)) shows that it is 
proportional to sin(0^ — 6o), and thus that our assumption that 6 = in the above model 
is not accurate. If we were to include both the inertial term and a driving term due to 
turbulence in the above model, then the equation would be equivalent to that of a damped 
driven pendulum, although we do not know how to represent the driving term. To better 
compare the data to the model, we calculated ^o(^O) ''') = [(^o{t + St/2) — 9o{t — 6t/2)]/St 
(where r represents a shift of the time axis), and binned the data based on 6o{t = 0) = 
[6'o((5t/2) + 9o{—5t/2)]/2. This yielded a rotation rate as a function of the delay time 
T after 9o had been reached. The above model neglects inertia, so we chose the delay 
time T = Tmax such that Oo{Oo,Tmax) is maximized to satisfy ()o{Tmax) = 0. Typically 
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Figure 31. (a): The coefficient of from (eo(6'o, Tmax)) = A{P) sin(6',3 - 6'o). Solid line: a fit 
of A{p) = A'\i3-/3o\ to the data for |/3| < 0.027. (b): The preferred orientation O/s obtained from 
the fit (open diamonds) and 9m from the peak of p(0o) (small circles). Data is for = 1.1 x lO'^" 
in the medium sample. 



Tmax ~ 100 s. The results for (^0(6*0, T„iaa;)) vs. 00 for i? = 1.1 X 10^° and /? = -0.0044 
rad are shown in Fig. 30. The error bars represent the probable error of the mean for 
each bin. A sinusoidal function (^0(^0, Tmax)) = Asin{9i} — 9q) was fit to the data to find 
the amplitude A{(}) of the deterministic rotation rate and preferred orientation 6/3(/3). 
It should bo noted that most of the time 6^ is near the preferred orientation (where 
{So{9Q,Tmax)) = and V{^()(0o, Tmaa;)) < 0), and thus the error bars are much smaller 
near this orientation in Fig. 30. For > 0.027, not all orientations are sampled and the 
fit to (^o(^Oi "^moa:)) is essentially a linear fit near Of). 

The fit parameter A{(3) is shown in Fig. 31a and 9fj{(}) (open diamonds) is shown in 
Fig. 31b, for several tilt angles /3 at i? = 1.1 x 10^° in the medium sample. Also plotted in 
Fig. 31b is the preferred orientation 6m (small circles) obtained from the peak of p(0o)- 
As we would expect, both methods of finding the preferred orientation agree with each 
other, so: 6m = ^/3, at least for tilt angles > 0.01 rad. Figure 31a shows a linear 
fit of A{(3) = A'\(3 - /3o| for < 0.026, which yields A' = 0.307 ± 0.002 rad/s and 
/3o = 0.00258 ± 0.00008 rad. These plots show that (^o(6'o, r^a^)) oc PsiniOp - 0o) as 
was predicted for f3 < 0.026, although the proportionality breaks down for larger /3. The 
calculated coefficient from the model of 0.40 rad/s is slightly larger than the fitted value 
of A', which we consider to be good agreement. This shows that the addded buoyancy 
in the tilted sample is directly responsible for the measured deterministic rotation rate 
for < 0.026. The fact that /3o 7^ is again due to the Earth's Coriolis force [Brown & 
Ahlers (2006)]. 

Finally we consider what effect the /?-dependence of the deterministic rotation rate has 
on reorientations. The minimum rotation rate for counting rcoricmtations for these data 
is 6min = 0.0254 rad/s, which is much larger than the deterministic forcing A = 0.0047 
rad/s found for (3 = —0.013, and it is even larger than the largest measured value for 
the deterministic forcing A = 0.015 rad/s. Thus wc conclude that the direct forcing 
from tilting the sample that is considered in the above model is too small to account 
for the 40% reduction in reorientations at (3 = —0.013 relative to the level sample and 
the complete suppression of reorientations within our resolution at larger values of 
There must be some other mechanism for the suppression of reorientations that we have 
not yet identified. 




9. Comparison with contemporary experiments 

In a set of contemporary experiments, Xi et al. (2006) report measurements of the 
orientation of the LSC for a RBC sample with T = 1, a w 5, and 10^ < i? < 10^°. 
To study short-term dynamics, they used particle image vclociiiicrtry to visiializc the 
horizontal fluid velocities near the top plate. From these measurements they calculated 
a spatially-averaged velocity and orientation for the LSC. For long-term measurements, 
thoy used a bead attached to a "fishing line" near the bottom plate to determine the 
orientation of the LSC. 

Many of the results reported by Xi et al. (2006) are similar to ours and provide an 
excellent complcnicnit to our measurements, as they use velocity mcasiirements and we use 
temperature measurements to quantify aspects of the LSC. Some of the results deserve 
special comment in relation to our work. 

Xi et al. (2006) report an oscillation of the LSC orientation aromid a preferred orienta- 
tion that they measured near the top and bottom plates. This oscillation was measured 
before by Funfschilling & Alilcrs (2004) on the basis of plume motion across the top and 
bottom plates. We find it as well, but only in the upper and lower rows of thermistors at 
heights 3L/4 and L/4. As seen by Funfschilling & Ahlers (2004), the oscillations in the 
top and bottom rows are out of phase with each other. This can be seen to some extent 
in Fig. 2 (we will report on it in detail elsewhere [Brown et al. (2006)]). We note that 
other aspects of the LSC azimuthal dynamics can depend quantitatively on the height 
at which they are measured, but so far this is the only process known to have such a 
qualitative height dependence. 

Xi et al. (2006) report that the LSC orientation remains locked near a preferred angle 
for the duration of an experimental run for a large majority of such runs. Curiously, this 
preferred orientation changes with each experimental run in a non-reproducible manner 
and thus cannot be due entirely to some geometrical asymmetry of their apparatus. It 
suggests a more complicated symmetry-breaking process. For comparison to our experi- 
ments we show two time scries of 6'o(i) over n.6 days at i? = LI x 10"'^^ in Fig. 32. The 
upper plot is for /3 = and the lower one is for (3 = —0.0044 rad. Consistent with the re- 
sults shown in Fig. 6, the /3 = data show a weak preference for some angle, presumably 
due to asymmetries such as the Coriolis force or very slight deformations of the side walls 
[Brown & Ahlers (2006)] ; but the preferred orientation is not nearly as severe as the one 
observed by Xi et al. (2006). The plot for /3 = —0.0044 rad looks qualitatively similar to 
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the time series reported by Xi et al. (2006), although the frequency of rotations through 
an entire revolution is much lower and p{0o) is still wider in our case. Since in our case the 
asymmetry is due to tilting the sample, we should not expect these statistical quantities 
to correspond exactly to those of Xi et al. (2006). 

Considering the differences between the two experiments, the cause for this preferred 
orientation may be the sapphire top plate in their experiment, with about 1/10 the con- 
ductivity of our copper plates. Direct numerical simulations (DNS) of the Boussinesq 
equations by Verzicco (2004) and experiments by Brown et al. (2005a) show that low- 
conductivity top and/or bottom plates reduce the heat transport across the sample. The 
DNS suggest that the low plate conductivity reduces the plume emission from the plates. 
Following Verzicco, we suggest that it is possible that the LSC put a significant thermal 
imprint on the sapphire plate. This thermal imprint would diffuse more quickly in a 
copper plate due to its higher conductivity. The thermal imprint in the sapphire plate 
presumably discourages cold plume emission, thus introducing a preferred orientation of 
the circulation corresponding to the initial orientation chosen by the LSC. If the thermal 
imprint remains longer than the duration of cessations and full azimuthal revolutions of 
the LSC, then the preferred orientation can remain so for the duration of the experi- 
mental run. This would also explain why the preferred orientation changes between runs, 
assuming the flow is stopped for some signiflcant length of time between them. 

To test this model, we examined data taken with aluminum top and bottom plates 
[Brown et al. (2005a)] in place of the usual copper ones. Aluminum has a thermal con- 
ductivity about 1/2 that of copper, which is still a factor of about four or five larger 
than that of sapphire, but since the aluminum plates slightly suppressed the overall heat 
transport relative to the copper-plate system, we hoped to see some influence also on the 
LSC orientation. In this case, the eight side-wall thermistors were glued to the outside 
of the side wall at a height L/2, so they were 0.31 cm away from the fluid. The thermal 
response time through the side wall is expected to be d'^/Kgw = 0(10) s. The locking 
of the orientation is compared in Fig. 33 by p{0o) for copper plates (open circles) and 
aluminum plates (solid circles). The distribution is more strongly peaked with aluminum 
plates and does not sample all angles as well as with copper plates , indicating that the 
orientation tends to be more strongly locked in place with aluminum plates. This sup- 
ports our model about the thermal imprint suppressing orientation changes of the LSC, 
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but a pair of experiments comparing copper and sapphire plates with a much larger ratio 
of conductivities would be more convincing. 

Xi et al. (2006) report finding a double-cessation event, in which there are two succes- 
sive cessations, and the second cessation returns the LSC to its original orientation. While 
we have not identified these events in our system, they are not necessarily inconsistent 
with our measurements. Since we count a cessation as occurring for the duration that 
the amplitude S remains below some threshold 5^, H & double cessation occurred in our 
system, we would most likely count it as a single cessation. Our measurements of 6 have 
a relative error of about 12% when S is near its average value, making it difficult for us 
to resolve double cessations. However, it is also quite possible that the double cessations 
are uniquely due to the suppression of plumes due to the thermal imprint on the sapphire 
plate, which could explain why the LSC is forced back into its original orientation. 

The average frequency of cessations (also counting double cessations) for the data 
reported by Xi et al. (2006) is loc = 1.7 day~^ with a probable error of the mean of 0.3 
day~^. This is consistent with our value of 1.5 ± 0.1 day~^. Considering the different 
types and locations of the measurements of the two experiments, and the sensitivity 
of cessation frequency to minor asymmetries as well as to how they are counted, it is 
somewhat fortuitous that the frequency of events in both experiments is in such good 
agreement. Xi et al. (2006) report an increasing frequency of cessations with increasing R, 
in contrast to our uniform frequency of cessations. This could be due to the criteria used 
for defining cessations: we used a minimum amplitude Si that changes with R, while Xi 
et al. (2006) use the velocity time series to find cessations. While we expect the velocity 
and temperature amplitudes to have similar behavior, we do not know if they should 
behave exactly the same during cessations, and thus we do not know if our methods of 
counting cessations are fully equivalent. 

The results from the analysis of crossings by Xi et al. (2006) (they referred to the 
events as "reversals") are qualitatively similar to ours. In particular, the shape of the 
probability distribution of the time intervals ti between crossings has the same shape 
(see our Fig. 17). Notably, the time interval at the crossover between the power-law de- 
pendence and the exponential dependence is about lOT in both cases. The characteristic 
decay time of the exponential region was significantly larger for Xi et al. (2006) (54T) 
than for us (32T). However, since this represents a typical long time interval between 
crossings, and the azimuthal motion was suppressed in their case, this difference is not 
surprising. 

10. Summary and conclusions 

We presented a broad range of measurements of the orientation of the LSC, includ- 
ing rotations and cessations. These events have not been well-studied experimentally or 
theoretically in the past, and we have very little physical understanding of how these 
phenomena occur. One important conclusion we can make is that when the LSC slows to 
a stop during a cessation, it loses information of its previous orientations and restarts at 
a random new orientation. We also found that both cessations and reorientations have a 
Poisson distribution in time, indicating the independence of successive events. We mea- 
sured the rate of change of the amplitude S during cessations. Its value can be compared 
with future dynamical theories of cessations. 

Sreenivasan et al. (2002) found that the time interval n between their successive 
events had a power-law probability-distribution with an exponent of minus one when ti 
was small, and was cut off exponentially at larger n . They, and in more detail Sreeni- 
vasan et al. (2004), interpreted the power law as indicative of self-organized criticality 
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(SOC) and attributed the exponential cut-off to a finite-size effect. In another paper Hwa 
et a/. (2005) proposed an analogy between the statistics of wind reversals on the one hand 
and that of fluctuations of the magnetization of the two-dimensional Ising model on the 
other. Again this analogy rested heavily upon the existence of a power-law probability- 
distribution for n. Our results for the statistics of the rotations and cessations that make 
up our reorientations are inconsistent with this interpretation and reveal a purely Pois- 
sonian probability distribution for our data which goes to a constant at small ti and 
drops off exponentially at large Ti; but when we include the relatively small-amplitude 
and high-frequency jitter of 9o{t) in the analysis by considering "crossings", then we 
reproduce the statistics observed by Sreenivasan et al. (2002). Thus we conclude that 
the SOC and the Ising-model analogy discussed by them and by Hwa et a/. (2005) can 
perhaps apply to the jitter of our measurements of 6*0 (t) but does not pertain to the 
reorientations {i.e. rotations and cessations) observed by us. We also have no reason to 
invoke a finite-size effect to explain a large- ti cutoff because the Poissonian statistics of 
our reorientations naturally yields exponential behavior at large ti . 

Both reorientations and cessations are found to be strongly suppressed in a tilted 
sample. The measurements could not be accounted for simply by the buoyancy of the 
boundary layer creating a preferred orientation of the LSC. This leaves an open problem 
for future work to answer. In addition it is unknown whether this reduction in the fre- 
quency of events applies to other types of asymmetries, such as a strong Coriolis force or 
more complicated geometries that arc common in geophysical systems. 

We presented data showing a statistical relationship between the rotation rate |^o| 
and the amplitude 5 of the LSC. This is interesting because of a phenomenon found in 
nature: reversals in the orientation of the Earth's magnetic field presumably as a result of 
convection reversal in the outer core are also known to be accompanied by a decrease in 
the amplitude of the resulting magnetic field [Glatzmaier et al. (1999)]. Previously Brown 
et al. (2005b) reported a correlation function of 5 and |^o| in which 5 tends to lead |^o| 
by about 6% of the turnover time. If this is extrapolated to the 400 year turnover time 
of the Earth's outer core, measurements of the Earth's magnetic-field amplitude could 
give several years' notice of a magnetic-field reversal. Convection in the Earth's core is in 
many ways different from our ideal convection experiment, so this extrapolation is highly 
speculative; but it is a subject worthy of future study if it could lead to a forecasting of 
events which would have a significant impact on human life. 
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12. Appendix: The velocity and temperature distributions in the 
plume layer 

The bimodal vertical velocity distribution near the side wall reported by Nicmela 
et al. (2001) led them to suggest that the LSC orientation switched between two opposite 
directions aligned with their sensors. We suggest that a bimodal local velocity distribution 
does not imply a bimodal probability distribution p{9o), and that the orientation could 
have been varying erratically through any azimuthal angle in their experiment. 

Recent particle velocimetry measurements in a cylindrical RBC sample of F = 1 and 
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Figure 34. (a): A model of the velocities in the outer shell of the LSC. (b): top view of the 
bottom plate, showing the length L of streamlines crossing the plate. Each streamline is referred 
to by the azimuthal angle 9 where it intersects with the side-wall boundary-layer. 

R = 7.0 X 10" provide some useful information [Sun et al. (2005)]. For the vertical velocity 
near the side wall in the plane aligned with the preferred orientation of the LSC (caused 
by a slight tilt), no reversals of the velocity direction were found. This agrees with our 
results for a tilted sample discussed in Sect. 8. In the plane orthogonal to this, many 
random velocity reversals were found, and they happened at the same time on opposite 
sides of the sample, so that there was always one side with up-flow and one side with 
down-flow Further, the magnitude of the velocity averaged over short time intervals (10 
min, without reversals) was about the same in both planes. This orthogonal-plane result 
is qualitatively similar to that of Niemela et al. (2001), although the latter did not have 
a tilted sample and did not know the preferred orientation of the LSC in their system. 
Since Sun et al. (2005) found no velocity reversals in the plane of the LSC, it follows that 
the observed reversals in the orthogonal plane were not likely to be global reversals of the 
LSC. Since the magnitude of the velocity was about the same in both planes, this suggest 
that the vertical velocity distribution near the side wall was nearly a step function with 
velocity +v at angles within ±7r/2 of and —v at all other angles, although realistically 
the velocity must smoothy transition from -\-v to —v over some small angular range. Since 
Niemela et al. (2001) only observed local velocity reversals, these could have mostly been 
small orientation changes in which the +/ — v interface moved past the sensors. If the 
sample of Niemela et al. (2001) was sufficiently symmetric, it is still likely that some of 
the measured events were rotations or cessations. 

The question remains whether this step- function velocity-distribution can be reconciled 
with our measurement of the shape of the temperature distribution around the side wall, 
which we found to be given by T{e) =To+6T with ST = S cos{0 - 9o) (see Fig. 1). Here 
we present a model of the thin shell between the boundary layer and the bulk, also called 
the plume layer. We assume that the average flow in this shell is parallel to the nearest 
wall of the container, and that the streamlines make complete loops and are parallel to 
the plane of the LSC as drawn in Fig. 34. This time average should be thought of as 
along enough to identify the LSC but short compared to the time between reorientations. 
Our calculation is based on the steady-state heat-transport equation 

(u- V)T = kV^T . 

We approximate the solution for the horizontal plume layers near the top and bottom 
plate of the sample, assuming that there the heat transfer to the plume layers is dom- 
inated by contact with the thermal boundary layers. The transport term is given by 
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{u- V)T = udT/dx « 2uST/Sx = 2u{9)6T{e)/\Dcos{0 - 0q)\ where we approximated 
the temperature gradient in the plane of the LSC by the horizontal temperature dif- 
ference 2ST{d) over the path length 6x = Dcos{9 — 9o). The heat flux to (from) the 
bottom (top) plume layer from the thermal boundary layer should not depend on 9, 
so V^T « +{—)constant. Putting these terms together yields the desired relationship 
between the horizontal velocity and temperature: 

u{9)5T{9) (X +(-) cos(6' - 9o) 

We found a horizontal temperature difference of the form 5T(9) = 5 cos{9 — 9o), which on 
the basis of our model implies that u{9) has a uniform magnitude but opposite direction 
near the top and bottom boundary layers. Assuming that viscous drag and other mech- 
anisms do not significantly affect the plume-layer velocity, following the streamlines in 
Fig. 34a implies that the vertical velocities along the side wall must have the same mag- 
nitude, resulting in the previously mentioned step-function velocity-distribution. This 
model ignores many aspects of the LSC, and thus cannot be a complete description of 
it. The model is meant to provide an understanding of why our measured temperature 
distribution around the side wall has a different shape than the velocity distribution 
suggested by the data of Sun et al. (2005). Additionally it justifies our suggestion of the 
step-function velocity- distribution in discussing the data of Niemela et al. (2001). With- 
out the qualitative features of the model described here it might be difficult to reconcile 
the results of the three experiments. 
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